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Abstract 


A  functional  representation  of  the  magnetic  field  contribution 
of  the  solar-wind-driven  electric  currents  on  the  magnetopause  of 
Jupiter  is  presented.  The  representation  accounts  for  the  precession 
of  the  magnetic  axis  and  the  consequential  diurnal  variation  of  the 
shape  of  the  magnetopause.  Existing  models  of  the  two  interior  field 
sources  (planetary  dipole  and  currents  of  trapped  plasma)  were 
incorporated  into  a  set  of  eight  three-dimensional  static  models. 
Each  model  assumed  a  10°  tilt  of  the  magnetic  axis  with  respect  to 
Jupiter's  rotation  axis.  This  set  includes  representations  published 
[Engle,  \992b]  for  a  tilt  toward  (a  =  0°)  and  away  from  (a  =  180°) 
the  Sun  in  the  noon-midnight  meridian  plane,  as  well  as  six 
complementary  analogous  models  for  azimuthal  angles  of  45°,  90°, 
135°,  225°,  270°,  and  315°.  The  representations  were  connected  by 
replacing  expansion  coefficients  with  continuous  functions  which 
replicated  each  of  the  eight  original  static  models  and  interpolated 
between  them.  The  result  is  a  time-dependent  functional 
representation  of  the  magnetic  potential  of  the  magnetopause  surface 
currents  in  terms  of  an  expansion  in  associated  Legendre  functions. 
This  model  yields  predicted  magnetic  field  components  for  any 
location  within  the  magnetosphere  as  Jupiter's  magnetic  axis 
precesses  during  the  course  of  a  Jovian  day  (9^  55^^  41^). 
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I.  Introduction 

In  December  1973,  on  board  the  Pioneer  10  spacecraft,  high- 
field  triaxial  fluxgate  magnetometers  provided  the  first 
instantaneous  vector  measurements  of  the  net  magnetospheric  field 
of  the  planet  Jupiter.  Without  further  processing  on  board  the 
spacecraft,  the  digitized  measurements  were  sent  directly  to  ground, 
where  raw  position  and  magnetic  field  data  were  translated  into  a 
Jupiter-centered  rectangular  coordinate  system.  Analysis  of  the  field 
data  proved  Jupiter's  magnetospheric  field  to  be  significantly 
different  than  had  been  anticipated  [Smith  et  al.,  1974].  An 
extrapolation  based  upon  a  detailed  knowledge  of  the  Earth's 
magnetosphere  was  consistent  with  the  intrinsic  planetary  field  of 
Jupiter.  Pioneer  and  later  Voyager  observations  suggested  a  dipole 
with  dipole  moment  of  about  420,000  nT-Rj^,  where  Rj  is  the  Jovian 

planetary  radius,  71433.6  km  (Figure  1). 


Figure  1 .  Dipole  approximation  of  the  intrinsic  planetary 

magnetic  field  of  Jupiter.  The  magnetic  axis  is  tilted  10° 
away  from  the  rotation  axis  toward  the  Sun. 


(0 


X  axis  is  toward  the  Sun  (in  units  of  Rj) 


However,  the  magnetosphere  was  much  larger  than  had  been 
predicted,  and  there  was  an  appreciable,  co-rotating  current  sheet 
located  in  and  near  the  magnetic  equatorial  plane.  This  current  sheet 
was  in  the  form  of  a  plasma,  a  gas  consisting  of  an  equal  number  of 
electrons  and  positive  charges.  This  current  sheet  extended  from 
about  5  Rj  (the  orbit  of  the  moon  lo)  to  nearly  100  Rj  at  the  time  of 

the  Pioneer  encounters.  Greater  solar  wind  pressure  and/or  reduced 
trapped-plasma  pressure  during  the  later  Voyager  spacecraft 
encounters  led  to  a  more  compressed  magnetosphere  and  an 
observed  current  sheet  extending  only  60  Rj  from  the  planet  [Ness  et 

al,,  1979].  In  each  case  the  current  sheet  increased  the  net  magnetic 
field.  The  intrinsic  planetary  dipole  and  the  currents  of  charged 
particles  trapped  near  the  planet  constitute  the  two  interior  field 
sources. 


Figure  2.  Jovian  model  magnetic  field  due  to  the  Voyager-observed 

interior  field  sources.  The  magnetic  axis  is  tilted  10°  from  the  planet's 
rotation  axis  toward  the  sun. 
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X  axis  is  toward  the  Sun  (in  units  of  Rj) 
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Figure  2  displays  the  magnetic  field  due  to  these  interior 
sources. 

The  outer  limit  of  Jupiter’s  magnetosphere  is  defined  by  the 
magnetopause,  a  literal  magnetic  field  boundary.  Only  a  weak 
interplanetary  magnetic  field  exists  outside  of  this  boundary. 
Located  where  the  incoming  solar  wind  is  deflected  by  the  interior 
magnetic  field,  the  magnetopause  is  sensitive  to  both  the  intensity  of 
the  solar  wind  and  to  the  specifics  of  the  interior  field  sources.  The 
solar  wind  is  a  stream  of  charged  particles,  mostly  protons  and 
electrons,  which  flows  from  the  sun.  This  high-speed  plasma, 
interacting  with  the  interior  field,  generates  currents  on  the 
magnetopause  which  enhance  the  field  intensity  interior  to  this 
boundary  but  cancel  the  field  outside  of  the  region.  On  the  planet’s 
solar  side,  the  magnetosphere  is  well-defined  by  the  magnetopause, 
but  in  the  anti-solar  direction  this  region  extends  millions  of 
kilometers  beyond  the  radius  of  the  planet. 

These  three  sources--the  intrinsic  planetary  dipole 
approximation,  the  co-rotating  current  sheet,  and  the  magnetopause 
surface  currents— describe  the  large-scale  magnetosphere  of  Jupiter. 
The  net  field  due  to  all  three  primary  field  sources  is  displayed  in 
Figure  3. 
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X-axis  is  toward  the  incident  solar  wind.  Distance  units  are  Jovian  radii 


II.  Preliminary  Models 

Magnetometer  measurements  from  four  spacecraft  flybys,  each 
along  single  line  trajectories  near  Jupiter's  orbital  plane,  have 
provided  the  foundation  for  several  global,  three-dimensional  static 
models  of  the  Jovian  magnetospheric  field.  There  is  no  unique  way 
to  model  the  observed  field.  Each  successive  model  attempts  to 
more  completely  and  more  accurately  represent  the  magnetosphere 
region  by  incorporating  the  effects  of  additional  physical  phenomena 
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and/or  the  conclusions  of  more  recent  and  reliable  investigations. 
Work  by  a  number  of  investigators,  including  Barish  and  Smith 
[1975],  Beard  and  Jackson  [1976],  and  references  cited  by  them,  led 
to  a  Jovian  magnetosphere  model  which  included  the  effects  of  the 
solar  wind  interaction  with  interior  sources  corresponding  to  Pioneer 
10  observations  [Engle  and  Beard,  1980].  Subsequent  Voyager 
observations  led  to  alternate  models.  Connerney  et  al.  recalculated 
the  planetary  dipole  [1982]  and  presented  a  co-rotating  current 
sheet  model  with  an  azimuthal  component  parallel  to  the  magnetic 
equator  [1981].  Engle  [1992a]  used  these  interior  sources  to  develop 
an  idealized  "no  tilt"  model  which  incorporated  the  effects  of  the 
solar  wind.  A  subsequent  work  [Engle,  \992b]  included  the  effects  of 
the  tilt  of  the  magnetic  axis  away  from  the  normal  to  the  incident 
solar  wind. 


III.  Modeling  the  Magnetopause  Surface  Currents 

The  incident  solar  wind  plasma  contains,  in  each  unit  of 
volume,  approximately  equal  numbers  of  positive  and  negative 
charges  traveling  with  the  same  average  velocity.  Each  particle, 
upon  crossing  the  magnetopause  and  entering  a  region  of  magnetic 
field,  experiences  a  magnetic  force  described  by 

F=qvxB  (1) 

where  q  is  the  electric  charge  on  the  particle,  v  is  the  particle’s 
velocity,  and  B  is  the  local  magnetic  field.  The  response  of  charged 
particles  to  this  force  is  a  deflection  of  positive  charges  in  one 
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direction  and  negative  charges  in  the  opposite  direction.  The  net 
effect  of  this  interaction  is  that  surface  current  elements  form  on  the 
magnetopause  boundary.  The  incident  neutral  particle  flux,  having 
no  net  current  associated  with  the  motion  of  its  particles,  has  the 
motion  of  its  constituents  reordered  in  such  a  way  as  to  produce  a 
ret  current  in  the  direction  of  the  deflection  of  the  positive  charges. 
The  magnitude  of  the  net  current  is  twice  the  magnitude  of  the 
current  associated  with  the  positive  charges  alone.  Associated  with 
every  current  is  a  related  magnetic  field,  as  described  by  the  Biot- 
Savart  Law, 

dB  =M£L  ^  ^  ^ 

^  ,  (2) 

where  Ids  is  a  small  current  element  and  r  is  the  distance  from  the 
current  element  to  the  field  point.  In  this  case,  the  configuration  of 
the  boundary  current  is  such  that  the  additional  magnetic  field 
contribution  acts  to  increase  the  field  inside  the  magnetopause  and  to 
cancel  the  magnetospheric  field  outside  the  boundary. 

Treating  the  many  incoming  particles  individually  is  an 
impractical  method  of  arriving  at  an  overall  surface  current 
configuration.  Instead,  these  equations  are  combined  with  concepts 
in  classical  thermodynamics  to  arrive  at  equations  of 
magnetohydrodynamics.  The  particular  principle  applied  to  calculate 
the  shape  of  the  magnetopause  is  to  balance  the  incoming  solar  wind 
flux  with  the  net  magnetic  scattering  effect  of  the  magnetospheric 
field.  The  magnetic  field  perpendicular  to  the  boundary  is 
presumably  negligible.  Equating  the  external  plasma  pressure  on  the 
boundary  of  the  magnetosphere  to  the  magnetic  pressure  of  the 
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internal  magnetic  field  at  this  location,  the  magnetohydrodynamic 
differential  equation  for  the  magnetosphere  boundary  is 

nxB  ±n  -9  =  0 

^  ^  (3) 

where  B  is  the  net  magnetic  field  at  the  boundary,  v  is  the  unit 
vector  in  the  direction  of  the  incident  solar  wind,  and  is  the  unit 
vector  normal  to  the  magnetopause  surface,  given  in  spherical  polar 
coordinates  by 

"  >*36  r  sin6  3(1)  (a\ 


with 


Kn  = 


1  +  (1|£)2  +  (_1_^)2 

*‘30  r  sin6  3<j) 


An  initial  magnetopause  surface  configuration  is  calculated  using 
equation  (3)  with  the  magnetic  field,  B,  represented  by  interior  field 
sources  alone.  After  an  initial  three-dimensional  zero-order  surface 
has  been  constructed,  the  next  step  is  to  calculate  the  currents  on 
that  surface  attributable  to  the  deflected  solar  wind  particles.  These 
zero-order  model  surface  current  elements  are  computed  using  the 
following  equation,  where  B  is  still  the  net  interior  sources'  magnetic 
field: 

J  =  n  X  B  =  ( — § - - - )r  -  (Bx  +  — ^ - )e  +  (Bg  +  — )0 


r  sin0  3<f)  36 


r  sin0 


The  magnetic  field  arising  from  the  computed  current  elements  is 
then  calculated  by  numerically  integrating  the  Biot-Savart  integral  of 


the  currents  over  the  entire  magnetopause  surface,  S  ,  as  follows: 

Bn,SC=/ 


(SsXB„.,)x 
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Subsequent  iterations  again  utilize  Equation  (3)  to  calculate  the 
magnetopause,  but  now  include  the  newly-calculated  contributions  of 
the  surface  currents  in  the  net  field,  B,  at  the  boundary.  This  cycle 
continues,  producing  higher-order  solutions,  until  no  significant 
change  in  the  magnetopause  shape  and  surface  current  elements  is 
observed  for  successive  iterations.  At  this  point,  convergence  is 
attained  and  the  model  is  termed  self-consistent. 

In  order  to  display  the  characteristics  of  the  net 

magnetospheric  model  it  is  necessary  to  represent  the  magnetic  field 
contributions  of  the  magnetopause  surface  currents  in  a  functionally 
continuous  manner  which  is  consistent  with  the  Biot-Savart  integrals 
over  the  surface  current  elements.  These  elements  are  generally 

computed  on  a  5°  by  5°  grid  of  the  magnetopause  surface.  To 

achieve  a  functionally  continuous  representation  for  the  entire  space 
inside  the  magnetosphere,  selections  of  calculated  magnetic  field 
components  for  a  large  number,  N,  of  representative  interior  grid 
locations  (r,0,<t))  are  fit  to  the  following  Schmidt-normalized  associated 
Legendre  function  expansion; 

*Viax  n 

Os(r,0,(|))  =  X  Z  ^nm  Pir(cos0)  cos(m(l)) 

n=l  m=0  (8) 

Bg(r,0,<j))  =  -  V  Og(r,0,(j)) 

<t>5  is  a  scalar  magnetic  potential  function  of  the  surface  current 

elements.  The  negative  gradient  of  the  potential  function  yields  the 
magnetic  field  due  to  this  source.  Thus,  given  a  set  of  interior  field 
values  and  grid  locations,  coefficients  Gnm  calculated  such  that 

the  sum  over  the  N  values  of  the  vector  residuals  squared  is 


minimized.  This  method  follows  the  general  technique  set  forth  by 
Mead  and  Beard  [1964]  in  their  modeling  of  the  Earth’s 

magnetopause. 

Self-consistent  models  for  both  the  Pioneer-observed  interior 
sources  and  the  Voyager-observed  interior  sources  have  been 
published  {Engle  and  Beard,  1980;  Engle,  1991;  Engle,  1992a, b].  In 
the  latter  Voyager  models,  a  self-consistent  magnetopause  and 
surface  current  calculation  was  fit  to  an  associated  Legendre 
functions  series  of  order  n  =  1  through  n  =  10.  This  expansion 
required  sixty-five  distinct  coefficients,  For  each  of  the 

magnetosphere  models  which  lay  the  groundwork  for  our 
investigation,  the  development  method  was  similar.  However,  until 
recently,  the  magnetopause  surface  calculation  for  each  configuration 
assumed  the  magnetosphere  surface  current  contribution  to  be  that 
achieved  by  the  idealized  "no-tilt"  calculation.  The  no-tilt 

assumption  greatly  expedited  the  achievement  of  self-consistency. 
Each  of  those  calculations  assumed  the  axis  of  all  the  interior  sources 
of  magnetic  field  to  be  perpendicular  to  the  incident  solar  wind. 

The  rotation  axis  of  Jupiter  is  tilted  3®  with  respect  to  the 
normal  to  the  planet's  orbital  plane,  and  the  magnetic  axis  makes  an 
angle  A,  =  10®  to  the  planet's  rotation  axis.  As  a  result,  the  magnetic 

axis  precesses  about  the  planetary  rotation  axis  during  each  diurnal 

rotation  of  period  9^  55®^  41^.  While  ignoring  the  small  tilt  of  the 
rotation  axis,  recent  model  calculations  have  incorporated  the  effects 
of  the  tilt  of  the  magnetic  axis  relative  to  the  incoming  solar  wind 
[Engle,  1992b].  A  precession  angle,  a,  is  a  useful  parameter  for 
describing  the  azimuthal  orientation  of  the  magnetic  axis  relative  to 


the  incident  solar  wind.  For  the  North  magnetic  pole  tilted  10° 
toward  the  sun,  a  =  0°;  tilted  away  from  the  sun  (half  of  a  rotation 
period  later),  a  =  180°.  Any  given  time  during  the  Jovian  day  may 
be  identified  with  a  corresponding  value  of  the  precession  angle,  a. 

Self-consistent  models  for  each  of  eight  values  of  a,  at  45° 
intervals,  have  been  constructed  by  means  of  the  iterative  method 
described  above.  Similar  processes  may  be  performed  for  any  static 
model  corresponding  to  any  particular  assigned  precession  angle 
orientation.  The  details  of  the  models  for  a  =  0°  and  a  =  180°  have 
already  been  published  [Engle,  1992b]. 

IV.  Developing  a  Time-dependent  Model 

The  limiting  feature  of  these  models  is  their  inherent  static 
nature.  Field  characteristics  at  a  given  location  can  be  derived  with 
considerable  accuracy,  but  the  procedure  demands  an  inordinate 
amount  of  time.  For  any  given  precession  angle,  the  Schmidt- 
normalized  spherical  harmonic  expansion  must  be  refit  to 
observations,  and  new  coefficients  of  the  harmonic  expansion  must 
be  calculated.  A  more  valuable  model  would  be  one  which  provides 
a  complete  model  description  of  the  magnetic  field  associated  with 
the  diurnally  varying  Jovian  magnetopause,  a  model  which  is 
functionally  continuous  in  time  as  well  as  in  position.  Such  a  model 
would  provide  a  magnetic  description  of  any  point  within  the 
magnetosphere  at  any  time  during  the  Jovian  day. 

A  simple  dipole  and  a  co-rotating  current  sheet  produce  field 
components  which  are  relatively  simple  to  model  in  a 
time-dependent  manner.  Both  sources  are  modeled  to  be  symmetric 


about  the  magnetic  axis.  In  accounting  for  the  tilt  of  this  axis,  the 
field  components  need  simply  to  be  translated  from  a  sinusoidally- 
varying  reference  frame  into  a  fixed  frame  defined  by  a  rotation  axis 
which  is  approximately  normal  to  the  orbital  plane  of  the  planet. 
The  crux  of  a  complete  time-dependent  model  lies  in  the  complexity 
of  coupling  the  iter:!tive  calculation  of  the  surface  current  field 
source  to  the  precession  of  the  magnetic  axis  about  the  rotation  axis. 

Associated  Legendre  functions  were  appropriately  used  to 
describe  the  magnetic  potential  due  to  currents  on  the 
magnetopause.  The  coefficients,  Gnm»  ^ct  to  weight  each  of  the 

associated  Legendre  functions  which,  collectively,  map  this  magnetic 
potential.  As  the  magnetic  axis  precesses,  these  functions  must  be 
combined  in  different  proportions  to  reflect  the  dynamics  of  the 
magnetopause  and  the  resulting  surface  currents.  These  coefficients 
introduce  the  variation  which  allows  this  field  to  change  with  the 
precession  of  the  magnetic  axis,  a. 

Since  each  complete  rotation  must  bring  the  field  back  to  its 
original  state  (assuming  the  constancy  of  parameters  such  as  solar 
wind  strength  and  sheet  current  density),  these  coefficients  assume  a 
functional  form  which  is  dependent  on  the  precession  angle  and  is 
necessarily  periodic  in  nature.  In  the  absence  of  well-developed 
evidence  to  the  contrary,  the  coefficients  are  also  assumed  to  be 
symmetric  for  precession  angles  which  are  symmetric  with  respect  to 
the  noon-midnight  meridian  plane. 

This  assumed  symmetry  was  a  fundamental  assumption  of  the 
Pioneer  and  Voyager  magnetometer  teams  and  underlies  the  work  of 
these  teams.  There  is  some  speculation  that  the  sheet  current 


undergoes  expansion  and  compression  as  it  moves  from  the  more 
closely-confined  day-side  to  the  more  expansive  night-side  of  the 
magnetosphere  and  back  again.  Other  related  considerations  are  that 
the  precession  of  the  current  sheet  may  lag  that  of  the  planetary 
dipole  by  some  fixed  azimuthal  angle  or  that  the  current  sheet  may 
undergo  differential  rotation  about  the  planet  (rotating  at  smaller 
angular  velocities  at  greater  distances  from  the  planet).  These  ideas, 
if  correct,  may  introduce  some  deviation  from  the  assumed 
symmetry  of  the  magnetopause  and  the  surface  current  contribution. 
Analysis  of  Ulysses  magnetometer  data  may  provide  further  insight 
into  these  phenomena  in  the  near  future. 

Included  in  the  complete  Voyager  Jovian  magnetosphere  model 
published  at  static  precession  angles  of  0®  and  180®  [Engle  \992b] 
are  the  corresponding  Legendre  function  expansion  coefficients  for 
both  orientations.  Using  the  magnetospheric  surfaces  and  associated 
grid  of  surface  current  elements  for  each  of  the  values  of  the 
precession  angle  a  =  0®,  45®,  90®,  135®,  180®,  225®,  270®,  and  315®, 
expansion  coefficients  were  computed  for  eight  orientations  of  the 
magnetic  axis  equally  spaced  over  the  course  of  a  Jovian  day.  In 
each  case,  Legendre  functions  of  order  one  through  ten  were  used  to 
maintain  the  self-consistency  of  the  surface  current  structure.  The 
result  was  eight  sets  of  sixty-five  values,  each  value  acting  to  weight 
a  particular  associated  Legendre  function  in  constructing  the 
magnetic  field  components  for  a  given  precession  angle.  This  data 
was  regrouped  into  sixty-five  sets  of  eight  values,  so  that  each  set 
now  represented  the  variation  of  one  particular  coefficient  over  the 
course  of  a  day.  Appendix  A  displays  these  coefficients. 
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Each  of  these  sets  was  then  graphed  as  G(n,m,a)  versus  a  and 
was  studied  for  an  inherent  pattern.  A  concise  mathematical 
representation  of  each  coefficient  would  enable  smooth  interpolation 
between  any  two  of  the  eight  adjacent  (in  a),  self-consistent 
magnetic  field  representations.  Most  of  these  self-consistent  model 
coefficients  were,  in  fact,  approximately  symmetrical  with  respect  to 
the  noon-midnight  plane.  One  such  example  is  displayed  in  Figure  4. 


Figure  4.  Variation  of  the  associated  Legendre  function  expansion 
coefficient  G(4,0,a)  as  a  function  of  the  precession  angle  of  the 
magnetic  axis. 
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In  some  cases,  the  patterns  were  more  difficult  to  represent 
functionally,  as  indicated  by  Figure  5. 

Figure  5.  Variation  of  the  associated  Legendre  function  expansion 
coefficient  G(5,5,a)  as  a  function  of  the  precession  angle 
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The  fitting  of  the  surface  current  representation  to  the  expansion 
coefficients  Gnm  fot  each  static  model  was  somewhat  sensitive  to  the 
particular  selection  of  computed  field  components.  For  this  reason, 
we  considered  each  data  point,  G(n,m,a),  to  have  some  range  of 
uncertainty.  For  each  graph,  a  unique  periodic  function  was 
developed  which  would,  ideally,  intercept  the  eight  self-consistent 
model  points.  However,  none  of  the  sixty-five  coefficients  was 
precisely  symmetric  over  the  course  of  a  day.  Often,  two  precession 
angle  orientations  which  were  symmetric  with  respect  to  the  noon- 
midnight  meridian  plane  yielded  slightly  different  values  for  a  given 
coefficient.  In  such  a  case,  the  corresponding  functional 

representation  assumed  an  average  value  of  the  coefficient  for  both 
orientations.  For  precession  angles  between  any  two  of  the  eight 
self-consistent  model  values,  each  function  traced  a  path  which 
would  reasonably  approximate  the  variation  of  the  coefficient  over 
the  corresponding  range  of  precession  angles. 

Initially,  different  functional  forms  were  used  to  model 
different  coefficients,  G(n,m,a).  However,  this  experimentation 
proved  a  fifth-order  cosine  series,  a  symmetric  and  periodic  function, 
to  be  capable  and  practical  for  describing  most  of  the  sixty-five 
expansion  coefficients.  Actually,  the  precise  form  of  the  function  was 
inconsequential.  The  aim  was  to  achieve  a  reasonable  mathematical 
description  of  each  coefficient,  and  clearly  there  was  no  unique 
solution.  The  final  functional  forms  of  each  coefficient  are  listed  in 
Appendix  B.  Using  these  time-dependent  coefficients,  the  Voyager- 
based  model  was  streamlined  and  the  laborious  procedure  detailed 
above  for  calculating  the  magnetopause  structure  for  a  given 


precession  angle  was  made  obsolete.  Fortran-based  computer 
programs  had  previously  been  developed  by  Engle  to  construct  the 
magnetopause  and  perform  the  least-squares  fitting  to  the  associated 
Legendre  functions.  With  the  functions  G(n,m,a)  substituted  into  the 
original  program  structure,  the  model  could  now  calculate  the 
Voyager  Jovian  magnetospheric  field  at  any  point  within  the 
magnetosphere  (r,0,<t))  and  at  any  precession  angle,  a.  This 
representation  adds  a  functional  estimate  of  the  magnetopause 
surface  current  contribution  to  the  two  interior  sources  of  the 
planet's  field. 


V.  Testing  the  Model 

The  new  model's  utility  in  describing  the  Jovian  magnetosphere 
environment  had  yet  to  be  demonstrated.  Clearly,  the  functions 
developed  to  model  the  coefficients  did  not  precisely  match  the  self- 
consistent  coefficients  point  for  point.  The  pertinent  assumption  was 
that  any  deviations  would  be  small  and,  averaged  over  the  sixty-five 
functions,  inconsequential.  Thus,  the  new  model  should  produce 
time-dependent  values  which  approximately  incorporate  and 
maintain  the  self-consistency  of  the  old  models. 

Utilizing  the  new  time-dependent  model  software,  a  second 
Fortran-based  program  generated  the  coordinates  of  individual 
magnetic  field  lines  emanating  from  the  surface  of  Jupiter.  To 
investigate  the  physical  nature  of  these  field  lines,  hundreds  of 
similar  lines  were  graphed,  originating  at  a  wide  spread  of  locations 
on  the  surface  of  Jupiter  for  numerous  values  of  the  precession  angle. 


Figure  6  shows  the  field  lines  ori  lating  in  the  noon-midnight 

meridian  plane  ((j)  =0°)  for  a  preces^  angle,  a  =  60°. 

Figure  6.  Noon-midnight  plane  projection  of  selected  magnetic  field  lines 
of  the  Jovian  magnetosphere  at  a  precession  angle,  a  =  60° 
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The  goal  was  to  ensure  that  families  of  magnetic  field  lines  predicted 
by  the  time-dependent  model  were,  in  fact,  physically  realistic.  A 
well-behaved  line  would  be  smooth  and  continuous  and  would  not 
cross  any  of  its  neighboring  field  lines.  Of  particular  interest  were 
those  lines  originating  near  the  cusp  regions.  The  two  cusp  regions 
(clearly  displayed  in  Figure  3)  are  the  "effective"  magnetic  poles  of 
the  magnetospheric  field.  Magnetic  field  lines  branch  out  from  the 
northern  cusp  and  return  to  the  southern  cusp.  The  true  magnetic 
poles  are  at  opposite  extremes  of  the  planet,  but  the  interaction  of 
the  solar  wind  causes  these  cusp  regions  to  be  swept  sunward  from 


the  poles.  Magnetic  field  lines  emanating  from  these  points  closely 
parallel  the  magnetopause  and,  if  modeled  incorrectly,  will  tend  to 
deflect  outside  of  the  magnetopause,  representing  nonphysical  field 
components  outside  of  the  magnetosphere  region. 

In  calculating  magnetic  field  lines,  those  which  extend  beyond 
70  Rj  from  the  planet  are  discontinued  at  r  =  70  Rj  to  reflect  only  the 

range  over  which  reasonably  reliable  field  predictions  are  expected. 
Voyager  2  observations  indicated  a  subsolar  point  of  about  60  Rj.  In 

the  tilted  models,  the  subsolar  point  is  the  radial  position,  measured 
in  the  magnetic  equatorial  plane,  which  is  closest  to  the  Sun.  This 
point  is  the  reference  position  for  the  spherical  harmonic  expansion 
of  the  Voyager-based  model.  The  magnetic  potential  within  the  day- 
side  magnetosphere  region  lends  itself  well  to  this  functional 
description,  and  the  night  side  is  modeled  reasonably  well  within 
about  60  -  70  Rj.  However,  as  the  night-side  magnetosphere  widens 

and  extends  millions  of  miles  in  the  anti-solar  direction  (and  hence, 
outside  of  a  spherical  boundary  of  r  =  60  Rj),  the  associated  Legendre 

function  representation  is  decreasingly  reliable  in  the  modeling  of 
the  magnetic  potential. 


VI.  The  Trajectory  Test 

Obvious  programming  errors  were  eliminated  under  the 
scrutiny  of  the  telltale  magnetic  field  line  traces.  The  qualitative 
result  was  a  physically  realistic  time-dependent  description  of 
Jupiter's  magnetospheric  field.  A  more  stringent  test,  of  course, 
would  address  the  quantitative  nature  of  the  model.  How  well  does 
it  predict  magnetic  field  values  in  the  vicinity  of  Jupiter? 


Unfortunately,  means  for  comparison  were  limited.  Until  the  recent 
Ulysses  mission,  only  four  spacecraft  had  traced  trajectories  through 
the  Jovian  magnetosphere,  each  equipped  with  magnetometers  which 
measured  the  magnetic  field  in  spacecraft-centered  rectangular 
coordinates.  Each  of  these  explorers  was  limited  to  a  pass  near  the 
orbital  plane  of  the  planet,  leaving  the  higher  latitudes  unexplored. 
Nonetheless,  Voyager  2  data  provided  a  reasonable  control  against 
which  to  check  the  accuracy  of  the  time-dependent  model.  This  data, 
in  particular,  was  chosen  because  the  model  assumed  the 
parameters  for  solar  wind  strength,  subsolar  point,  and  sheet  current 
density  which  were  consistent  with  Voyager  2  observations. 

On  board  Voyager  2,  field  component  measurements  were 
taken  every  sixty  milliseconds.  For  practical  purposes,  these 
measurements  were  averaged  over  forty-eight  second  intervals, 
1800  data  points  per  Earth  calendar  day.  Position  determinations 
were  made  separately  and  followed  roughly  one  second  after  the 
corresponding  averaged  field  readings.  This  data  was  translated  into 
Jupiter-centered  rectangular  coordinates,  with  the  x-axis  in  the 
direction  of  the  Sun  and  the  z-axis  perpendicular  to  the  planet's 
orbital  plane.  Coupling  the  lagging  position  measurement  to  the 
leading  field  measurement,  these  readings  were  treated  as  though 
they  had  been  taken  simultaneously.  This  approximation  did  not 
introduce  any  significant  error. 

Magnetometer  readings  on  5  and  6  July,  1979,  revealed  that 
the  spacecraft  was  experiencing  Jupiter’s  bow  shock,  evidenced  by 
magnetic  "noise"  on  the  order  of  5  -  10  nT.  The  bow  shock  is  a 
supersonic  shock  wave,  a  dynamic  buffer  region  which  envelops  the 
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deflected  solar  wind  particles  and  separates  the  interspace  medium 
from  Jupiter's  magnetosphere.  The  bow  shock  is  characterized  by 
turbulent  flow  of  the  deflected  solar  wind  particles.  The  spacecraft 
crossed  the  magnetopause  early  on  the  6  July.  Unfortunately,  data 
from  the  Voyager  on  7  July  never  reached  Earthbound  receivers. 
This  limits  our  control  data  to  a  range  between  8  July,  when  Voyager 
2  was  inbound  at  about  33  Rj  from  Jupiter,  and  14  July,  when  the 
Voyager  was  outbound  on  the  anti-solar  side  at  about  72  Rj  from 

Jupiter.  Figure  7  shows  the  trajectory  of  Voyager  2  as  the  position 
data  was  received  by  earthbound  receivers. 


Figure  7.  1979  Jupiter-encounter  trajectory  of  Voyager  2  in  planetocentric 
orbital  coordinates.  The  x-y  plane  lies  in  the  orbital  plane  of  the  planet. 
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X  axis  is  toward  the  Sun  (in  units  of  Rj ) 


A  few  other  clearly  erroneous  position  and  field  measurements  were 
excluded  from  this  display  and  from  the  following  data  analysis.  In 
addition,  reports  of  initial  Voyager  2  observations  [Ness  et  al,  1979] 
indicate  that  some  data  from  8  July  may  be  tainted  by  an  apparent 
tumbling  of  the  spacecraft  which  was  unaccounted  for  in  subsequent 


reference  frame  translations.  The  magnitude  of  the  net  field 
observed  on  this  day  is,  of  course,  still  reliable,  but  components  of 
the  field  may  not  be  accurately  decomposed. 

The  reported  position  data  were  translated  into  Jupiter- 
centered  spherical  coordinates  and  the  resulting  data  set  was  used  as 
position  input  to  the  time-dependent  magnetosphere  model.  Again, 

Fortran-based  programs  were  used  to  handle  the  enormous  volume 
of  data.  The  model  calculated  the  magnetic  field  components  (Bj.,  B0, 

B^jj)  in  spherical  coordinates  for  each  data  point  along  the  Voyager 

trajectory.  Incorporated  in  the  model  is  the  rotational  rate  of  Jupiter, 
360°/(9^  55*”  41*),  which  equates  to  0.483480®  every  forty-eight 
seconds  (each  position  interval). 

It  was  necessary  to  assume  an  initial  precession  angle,  a 
reference  orientation  for  the  first  trajectory  position  on  8  July.  At 
each  subsequent  trajectory  point,  the  model  would  update  the 
precession  angle  based  on  the  planetary  rotation  rate  and  would  then 
calculate  the  magnetic  field  components.  An  initial  precession  angle 
ai  =  210®  at  00^00*”  8  July  was  chosen  such  that  the  magnetic  axis 

was  tilted  directly  toward  the  sun  at  00^00*”  of  the  day  of  Jupiter's 
closest  approach  to  the  planet,  9  July.  This  arbitrary  choice  was  a 
matter  of  convenience.  The  ensuing  field  calculations  provided 
grounds  for  calibrating  this  assumed  initial  precession  angle  with 
spacecraft  observations. 

Next,  the  Voyager-observed  field  measurements  were 
translated  into  a  form  which  was  comparable  with  model 
calculations.  A  Fortran-based  program  utilized  the  same  spacecraft 
position  data  file  which  was  employed  by  the  time-dependent  model. 
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but  it  also  converted  the  corresponding  observed  field  component 
measurements  (Bj.,  B0,  B^j^)  into  a  Jupiter-centered  spherical 

coordinate  system  with  <j>  =  0®  toward  the  Sun. 

Because  the  Voyager  2  trajectory  paralleled  Jupiter's  orbital 
plane,  the  spacecraft  encountered  the  planet  at  locations  in  and  near 
the  magnetic  equatorial  plane.  The  precession  of  the  magnetic  axis 
caused  the  co-rotating  current  sheet  to  repeatedly  cross  the 
spacecraft  trajectory.  Hence,  the  Voyager’s  Jupiter-centered 
coordinates  alternated  between  the  northern  magnetic  hemisphere 
and  the  southern  magnetic  hemisphere  during  the  course  of  each 
day.  Referring  to  the  field  line  plot  in  Figure  3,  it  is  clear  that 
magnetometer  readings  of  the  component  Bj.  should  indicate 

alternating  positive  and  negative  values  as  the  Voyager  2  entered 
the  northern  and  southern  magnetic  hemispheres,  respectively. 
When  the  magnetic  axis  was  tilted  away  from  Voyager  2,  la-<(>l  >  90®, 
field  lines  were  directed  radially  inward  to  the  planet,  while  a  tilt 
toward  the  Voyager,  la-<l)l  <  90®’  caused  f  eld  lines  to  point  radially 
outward. 

The  nearly  periodic  nature  of  the  Bj.  component  provided  a 

simple  means  to  calibrate  an  initial  precession  angle.  The 
assumption  that  aj  =  210®  had  resulted  in  Bj.  component  calculations 

which  were  almost  exactly  180®  out-of-phase  with  spacecraft 
observations.  Rotating  the  magnetic  axis  forward  180.2®  in  its 
precession  put  model  calculations  of  the  current  sheet  crossings  in 
phase  with  spacecraft  observations.  Setting  ai  =  30.2®  at  00^  OO”’  on 

8  July,  the  model  was  again  used  to  predict  the  field  components 
along  the  following  seven  days'  journey  through  the  magnetosphere. 
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The  model  values  and  observed  values  of  the  Bj.  component  on  9  July 
are  overlaid  in  Figure  8. 


Figure  8.  Model  and  observed  magnetic  field  components  directed 
radialiy  outward  from  Jupiter  on  9  July,  1979. 


Time  (in  units  of  48  seconds) 


The  agreement  between  theory  and  experiment  is  remarkable. 

The  model  accurately  portrays  the  sheet  current  crossings  and  also 
predicts  values  for  Bq  which  nearly  coincide  with  the  observed 

values.  While  the  variation  in  Bj.  closely  matches  observation,  the 

model  does  predict  a  greater  range  of  magnitudes  than  are  observed 
for  this  field  component.  Curiously,  B^j,  seemed  to  be  180°  out  of 

phase  with  observations,  although  the  much  smaller  magnitude  of 
the  B^  component  made  comparisons  more  difficult  and  less 


meaningful. 


VII.  Calibration  and  Some  Deductions 

Voyager  2  observations  were  useful  in  substantiating  what  was 
now  the  first  mathematical  description  of  the  Jovian  magnetosphere 
to  incorporate  the  continuous  diurnal  variation  of  the  magnetopause 
current  contributions.  Model  calculations  proved  to  be  reasonably 
consistent  with  spacecraft  observations.  Through  a  closer 
comparison  of  these  data  sets,  an  assessment  of  various  model 
parameters  can  be  made,  leading  to  improvements  in  the  model  and 
deductions  about  the  actual  magnetic  field  source  structure.  Several 
attempts  were  made  to  optimize  the  fit  of  the  model  to  Voyager 
observations  through  the  variation  of  model  parameters.  Some 
conclusions  derived  from  these  efforts  follow. 

Figure  9  displays  the  model  values  and  the  observed  values  of 
the  Bj.  component  over  the  entire  trajectory  data  set,  8-14  July. 


Figure  9.  (a)  Model  magnetic  field  components  directed  radially  outward 
from  Jupiter  from  8-14  July,  1979. 
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Figure  9  (b)  Voyager-observed  magnetic  field  components  directed  radially 
outward  from  Jupiter  from  8-14  July. 


The  precise  agreement  of  the  model's  predicted  current  sheet 
crossings  with  the  Voyager-observed  crossings  (where  Bj.  =  0,  we 

infer  that  Voyager  2  crossed  the  center  of  the  current  sheet) 
substantiates  the  concept  of  a  rigid  current  sheet.  This  is  significant 
because  there  has  been  some  speculation  in  recent  years  that  this 
current  sheet  is  not  rigid  but  actually  experiences  differential 
rotation  about  the  magnetic  axis.  This  would  entail  the  current  sheet 
rotating  with  slower  angular  velocity  as  distance  from  the  planet 
increases.  If  this  were  the  case.  Voyager  2  would  not  have  recorded 
the  consistent  periodicity  of  the  sheet  crossings  which  is  indicated  by 
the  both  the  spacecraft  data  and  the  rigid  current  sheet  model 
(incorporated  in  the  time-dependent  model).  The  fact  that  it  did 
measure  this  periodicity  is  confirmation  of  the  rigidity  of  the  sheet 


current  structure. 
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On  12  July,  the  Voyager  2  made  the  subtle  transition  into  the 
magnetotail  region  of  Jupiter,  where  it  remained  for  the  following 
two  days  of  our  data  set.  Stretching  millions  of  kilometers  in  the 
anti-solar  direction,  the  magnetotail  region  is  defined  by 
electromagnetic  phenomena  which  are  not  properly  addressed  by  the 
magnetosphere  model.  Within  a  radial  distance  of  about  35  Rj  in  the 

anti-solar  direction,  current  sheet  crossings  occurred  nearly 
coincident  with  the  spacecraft  traversal  of  the  magnetic  equatorial 
plane.  As  the  outbound  spacecraft  approached  the  pre-dawn 
magnetosphere,  however,  south-to-north  sheet  crossings  occurred 
with  less  regularity.  This  phenomenon  has  been  previously  noted 
and  is  generally  understood  in  terms  of  a  "periodic  rocking  of  the  tail 
current  sheet  about  the  longitudinal  axis  of  the  tail,  as  Jupiter 
rotates".  [Ness  et  al,  1979]  The  magnetotail  region  is,  in  fact,  a 
distinctly  different  region  of  magnetic  activity  which  is  not 
addressed  by  models  of  the  Jovian  magnetosphere.  In  this  region, 
larger  discrepancies  were  noted  between  model  calculations  and 
Voyager  2  observations  than  during  previous  days  when  the 
spacecraft  was  nearer  the  planet. 

In  an  attempt  to  optimize  the  fit  of  the  model  to  Voyager 

observations,  maximum  likelihood  theory  was  applied  to  the 

variation  of  several  model  parameters.  Variation  of  the  initial 

precession  angle  led  to  further  deductions  about  the  Jovian  magnetic 

field  source  structure.  Based  on  the  visually  calibrated  initial 

2 

precession  angle,  aj  =  30.2®,  X  calculations  of  B^,  B0,  B^  and  B^q^^I 
were  performed  using  methods  consistent  with  those  outlined  bv 
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Taylor  [1992].  values  for  each  field  component  5C^0,  X^^) 

were  calculated  for  each  of  several  assumed  initial  precession  angles. 

Taken  individually,  these  calculations  meant  very  little. 
However,  a  comparison  of  the  X  values  from  each  of  several  closely- 
spaced  initial  precession  angles  provided  a  relative  comparison  of  the 
model’s  ability  to  match  the  observed  results  for  a  given  aj.  This 

could  make  possible  a  more  accurate  calibration  of  the  initial 
precession  angle  by  seeking  a  minimum  value  of  X  over  the  360® 

range  of  values  for  the  initial  precession  angle.  The  results  were 
surprising.  X  r  minimized  for  aj  =  67.2  ±  0.5®.  This  angle  is 

displaced  37®  from  the  visually  calibrated  initial  precession  angle,  aj 

=  30.2® .  A  careful  inspection  of  the  software  revealed  no 

programming  errors,  but  neither  can  this  discrepancy  be  positively 

attributed  to  any  other  cause. 

2 

The  minimum  X  9  value  occurred  for  an  angle,  aj  =28.2  ±  0.5®, 

2 

39®  displaced  from  that  for  the  minimum  X  r  calculation.  Near 

Jupiter's  orbital  plane,  the  69  component  is  strongly  influenced  by 

the  orientation  of  the  magnetic  axis.  The  Bj.  component,  on  the  other 

hand,  is  overwhelmingly  the  product  of  the  co-rotating  current  sheet. 

2 

Disparate  minimum  X  calculations  could  indicate  that  the  two 

internal  field  sources— the  planetary  dipole  and  the  co-rotating 

current  sheet— are  out  of  phase.  In  other  words,  the  current  sheet 

might  lag  the  magnetic  axis  in  its  precession. 

2 

Inferences  drawn  from  the  X  calculations  are  by  no  means 

conclusive.  Reservations  about  these  deductions  stem  from  the 
observation  that  the  minimum  X  r  value  does  not  correspond  to  the 
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visually  calibrated  initial  precession  angle.  Nonetheless,  the  ideas 
which  are  presented  demand  consideration  as  genuine  physical 
phenomena  of  the  Jovian  magnetosphere  system.  The  time- 
dependent  nature  of  this  model  makes  it  particularly  useful  in  the 
analysis  of  these  and  similar  experimental  observations. 

VIII.  Future  Tests 

While  a  new  magnetosphere  model  may  fit  observed  data 
somewhat  better,  this  is  no  guarantee  that  the  model  itself  is  better. 
The  essential  question  of  a  model's  utility  is  how  well  it  predicts  field 
properties  at  points  not  on  the  Pioneer  and  Voyager  trajectories, 
nearer  Jupiter  or  at  higher  or  lower  latitudes,  for  example. 

Average  sheet  current  density,  subsolar  point  and  other  model 
parameters  were  reported  for  the  Pioneer  10  and  11  flybys  as  well 
as  both  Voyager  flybys.  Substituting  Pioneer-deduced  parameters 
for  those  of  the  Voyager  mission  which  currently  define  the  time- 
dependent  model,  we  should  achieve  a  reasonable  approximation  of 
the  magnetic  environment  encountered  in  1973  and  1974  by  Pioneer 
10  and  11. 

More  recently,  in  1992,  the  spacecraft  Ulysses  encountered 
Jupiter  en  route  to  its  sweep  of  the  polar  regions  of  the  Sun.  The  data 
collected  at  Jupiter,  which  has  not  yet  been  released  for  public  use, 
will  provide  another  control  against  which  to  gauge  the  value  of  the 
Voyager-based  time-dependent  model.  Because  Jupiter  was  used  in 
a  gravity  boost  maneuver  to  deflect  Ulysses  out  of  the  solar  ecliptic 
plane  and  into  a  polar  orbit  of  the  Sun,  measurements  from  Ulysses 
spanned  a  much  wider  range  of  Jovian  latitudes  and  provided  a  more 
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complete  sampling  of  the  planetary  magnetic  field  than  any  previous 
mission  to  the  planet.  Data  from  this  mission  can  be  gainfully 
compared  to  the  predictions  of  our  time-dependent  model  along  the 
spacecraft  trajectory. 


IX.  Application  of  the  Model 

Perhaps  the  most  interesting  application  of  the  new  time- 
dependent  model  involves  the  upcoming  July  20,  1994  collision  of 
comet  Shoemaker-Levy  9  with  Jupiter.  A  team  headed  by  Principal 
Investigator  Renee  Prange,  of  the  Institut  d'Astrophysique  Spatiale 
at  the  University  of  Paris,  ORSAY,  France,  has  been  awarded 
telescope  time  on  the  Hubble  Space  Telescope  during  the  time  of  the 
comet's  pass  through  the  Jovian  magnetosphere.  This  team  will 
image  the  enhanced  particle  activity  which  results  when  ions  from 
the  comet's  corona  encounter  the  magnetic  field  of  Jupiter. 

As  the  comet  crosses  the  Jovian  magnetopause,  the  charged 
particles  which  constitute  the  comet’s  corona  will  be  subject  to  the 
same  magnetic  force  which  describes  the  solar  wind  deflection  at  the 
magnetopause  boundary.  These  particles  will  undergo  a  cyclotron 
process,  spiraling  up  and  down  the  magnetic  field  lines  which  cross 
the  comet  ephemeris,  or  trajectory. 

At  the  request  of  Prange,  and  based  on  the  predicted  comet 
trajectory  positions  which  she  provided,  Engle  originally  used  the 
static  self-consistent  model  to  calculate  the  magnetic  field  lines 
anticipated  to  lie  along  the  comet  ephemeris.  More  recently,  the 
time-dependent  model  has  been  applied  to  the  trajectory  for  a  wide 
range  of  magnetic  axis  precession  angles  to  locate  and  map  these 


predicted  magnetic  field  lines.  The  results  have  been  sent  to  the 
observation  team. 

Figure  10  displays  two  characteristic  magnetic  field  line 
calculations  corresponding  to  the  17.5  July  (12^  00”^  on  17  July)  and 
18.2  July  (04^  48^  on  18  July)  positions  of  Shoemaker-Levy  9. 


Figure  10.  Comet  Shoemaker-Levy  9  ephemeris  with  selected 

^  calculated  magnetic  field  lines  for  a  =  a. 

oc 
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Y-axis  is  toward  the  dusk  (in  units  of  Rj) 


The  first  line  encountered  on  the  inbound  trajectory  fails  to  close  on 
the  planet,  indicating  that  the  comet  is  close  to  crossing  the 
magnetopause.  The  magnetosphere  model  is  only  valid  for  points 
within  the  magnetopause.  For  points  just  outside  the  boundary,  the 
model  breaks  down  and  predicts  field  lines  which  fail  to  close.  The 
second  line  does  eventually  close  on  the  planet,  indicating  that,  at 
this  location  and  precession  angle,  the  comet  will  have  entered  the 
magnetopause.  These  and  similar  field  predictions  should  help  the 


team  to  anticipate  the  times  and  regions  of  enhanced  particle  activity 
within  the  magnetosphere. 

The  expected  particle  activity  is  of  particular  interest  to 
physicists  studying  the  comet's  composition.  Current  theory  places 
the  origin  of  comets  in  a  spherical  shell,  or  cloud,  some  50,000  AU,  or 
7.48  (10^2)  km,  from  the  Sun  [Chaisson  and  McMillan,  1993]. 
Because  this  cloud  is  believed  to  be  left  over  from  material  that 
condensed  to  form  the  solar  system,  the  composition  of  comets 
should  reflect  the  composition  of  the  original  solar  nebula.  Imaging 
the  evaporated,  ionized  constituents  of  the  comet  Shoemaker-Levy  9 
should  yield  significant  information  about  this  composition. 

X.  Conclusions 

Among  the  benefits  of  the  time-dependent  representation  is  its 
ability  to  produce  a  "living"  model  of  the  complete  magnetosphere  of 
Jupiter,  including  those  regions  which  are  near  to  the  magnetopause 
and  those  which  vary  substantially  in  their  field  properties  with  the 
precession  of  the  magnetic  axis.  The  model  stands  up  to  the 
qualitative  requirement  that  it  produce  physically  realistic  magnetic 
field  lines.  It  also  yields  field  component  calculations  which  are 
quantitatively  consistent  with  Voyager  2  magnetometer  readings. 

Clearly,  the  time-dependent  magnetosphere  model  is  already 
providing  a  service  to  the  space-science  community.  The  model  has 
been  used  to  provide  insight  into  the  upcoming  interaction  of  comet 
Shoemaker-Levy  9  with  the  Jovian  magnetosphere.  In  upcoming 
months,  this  model  may  also  serve  an  important  role  in  the  analysis 
of  the  recent  Ulysses  spacecraft  encounter  data  and  of  the  upcoming 
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Galileo  observations.  In  these  and  other  future  encounters  with 
Jupiter,  a  time-dependent  model  can  render  far  less  formidable  the 
task  of  reconciling  large  volumes  of  data  in  the  form  of  magnetic 
field  component  measurements  to  models  which  include  all  the  basic 
salient  features  of  the  magnetospheric  environment.  This  allows  a 
more  reliable  separation  of  phenomena  between  those  which  are 
inherently  characteristic  of  a  dynamically  stable,  albeit  diurnally 
varying,  magnetosphere  and  those  which  are  transient  perturbations 
in  the  magnetospheric  field. 

Developed  as  an  extension  of  existing  Voyager-based 
representations,  the  time-dependent  model  reflects  the  effects  of 
nearly  maximum  solar  wind  pressure  conditions.  However,  the 
model  should  provide  a  guide  to  analogous  representations  for  the 
magnetopause  surface  current  contributions  resulting  from  minimum 
solar  wind  pressure  conditions  (as  observed  by  the  Pioneer 
spacecraft)  and  solar  wind  pressure  conditions  between  the  two 
extremes. 

A  time-dependent  representation  should  be  valuable  to 
experimenters  who  directly  measure  characteristics  of  the  plasma 
particles  encountered  during  spacecraft  flybys.  These  particles  are 
charged  and  are,  thus,  sensitive  to  the  magnetic  field  environment 
through  which  they  move.  Knowledge  of  the  magnetic  field 
environment  experienced  by  those  particles  just  prior  to  their 
observation  is  essential  to  being  able  to  draw  inferences  regarding 
the  overall  population  of  the  plasma  environment  from  the  small 
sample  actually  observed. 
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Appendix  A 


Associated  Legendre  Function  Expansion  Coefficients  G(n,m)  for  Each 
of  Eight  Distinct  Values  of  the  Precession  Angle,  a 


a=0® 

a=45® 

a=90® 

a=135® 

G(  1.  0) 

-0.58087701 

-0.58464450 

-0.59458369 

-0.61738110 

G(  1,  1) 

-0.08373246 

-0.10518049 

0.01739297 

0.09397754 

G(  2,  0) 

0.01284187 

-0.00414274 

-0.00586055 

0.01140946 

G(  2,  1) 

-0.34831586 

-0.31549028 

-0.31097990 

-0.34939179 

G(  2,  2) 

-0.04608620 

-0.02888788 

0.00147852 

0.04055239 

G(  3,  0) 

-0.07753199 

-0.09483284 

-0.06431898 

-0.08158017 

G(  3,  1) 

-0.04682504 

-0.04671578 

-0.00771004 

0.05568385 

G(  3,  2) 

-0.07688551 

-0.05731300 

-0.05547160 

-0.07043299 

G(  3,  3) 

-0.01007293 

-0.00457890 

-0.00789330 

0.01073236 

G(  4.  0) 

0.04050628 

0.02233014 

-0.00737017 

-0.02311350 

G(  4,  1) 

-0.14027660 

-0.14176409 

-0.11574018 

-0.13151264 

G(  4,  2) 

-0.02936912 

-0.04489575 

0.00505519 

0.04719570 

G(  4,  3) 

-0.00192175 

-0.00283145 

0.00963413 

-0.00178103 

G(  4,  4) 

-0.00026736 

-0.00609120 

0.00217616 

0.00336623 

G(  5,  0) 

-0.00440191 

-0.00285284 

-0.00698109 

-0.00666580 

G(  5,  1) 

-0.01023021 

-0.02333847 

-0.00371074 

0.01137838 

G(  5,  2) 

-0.06185311 

-0.07630672 

-0.07628365 

-0.05615139 

G(  5,  3) 

0.00137553 

-0.01738635 

-0.00019316 

0.03860091 

G(  5,  4) 

-0.00061794 

-0.01071075 

-0.00189747 

-0.00169872 

G(  5,  5) 

0.00102967 

-0.00419982 

0.00550852 

0.00543485 

G(  6,  0) 

0.02902861 

0.02742495 

0.00004615 

-0.01623161 

G(  6,  1) 

-0.06302944 

-0.04277107 

-0.04756158 

-0.04363710 

G(  6,  2) 

-0.02455459 

-0.02079814 

-0.00113889 

0.02152534 

G(  6,  3) 

-0.00076756 

-0.02664453 

-0.03697173 

-0.02952027 

G(  6,  4) 

-0.00107624 

-0.00697785 

-0.00109286 

0.01247352 

G(  6,  5) 

-0.00561647 

-0.00695962 

-0.01348485 

-0.01886771 

G(  6,  6) 

0.00107300 

-0.00135418 

0.00238792 

0.00654377 

G(  7,  0) 

0.02305909 

-0.00585131 

0.01376711 

-0.00308078 

G(  7,  1) 

0.01690746 

0.00244996 

-0.00784696 

-0.01351918 

G(  7,  2) 

-0.04826577 

-0.02443572 

-0.03400584 

-0.02146269 

G(  7,  3) 

-0.00287149 

-0.01368837 

-0.00230085 

0.03535789 

36 
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G(  7,  7) 

-0.00210364 

0.00251089 

0.00473227 

-0.00330430 

G(  8,  0) 

-0.02000607 

-0.00113438 

0.00118011 

0.00501335 

G(  8,  1) 

-0.00000762 

0.00000660 

-0.00000005 

-0.00008275 

G(  8,  2) 

0.00019176 

-0.00007228 

0.00000001 

0.00002353 

G(  8,  3) 

-0.02702314 

-0.01217487 

-0.02555243 

-0.01238130 

G(  8,  4) 

-0.00447119 

0.00616529 

0.00269019 

-0.00550030 

G(  8,  5) 

0.00695114 

0.00245640 

-0.00477784 

-0.00017510 

G(  8,  6) 

-0.00974452 

0.00796753 

-0.00120304 

0.00479962 

G(  8,  7) 

0.00226276 

-0.00239843 

-0.00091922 

0.00006179 

G(  8,  8) 

-0.00202971 

0.00056569 

0.00433292 

-0.00297690 

G(  9,  0) 

-0.01033668 

-0.01000137 

-0.00067857 

-0.03044157 

G(  9,  1) 

0.00152764 

-0.00010562 

-0.00048214 

-0.00004036 

G(  9,  2) 

-0.00004681 

0.00005837 

-0.00021124 

0.00007267 

G(  9,  3) 

-0.00476104 

0.00161330 

0.00228422 

-0.00166715 

G(  9,  4) 

-0.00464545 

-0.00275073 

-0.00940238 

-0.00413898 

G(  9,  5) 

-0.00280360 

0.00289581 

-0.00217455 

-0.00216732 

G(  9,  6) 

-0.00207030 

0.00070886 

-0.00249742 

-0.00017582 

G(  9,  7) 

-0.00409556 

0.00299261 

-0.00062696 

0.00255423 

G(  9,  8) 

0.00127051 

-0.00028279 

0.00041439 

-0.00059604 

G(  9,  9) 

-0.00326805 

0.00052126 

0.00333244 

-0.00092754 

G(10,  0) 

-0.00094684 

0.00003501 

0.00017572 

0.00002547 

G(10,  1) 

-0.00025991 

-0.00016511 

-0.00009792 

-0.00031670 

G(10,  2) 

0.00112341 

-0.00080817 

0.00122698 

0.00013672 

G(10,  3) 

-0.00230579 

-0.00099702 

-0.00252340 

-0.00023840 

G(10,  4) 

-0.00318032 

-0.00104373 

0.00281916 

-0.00048342 

G(10,  5) 

0.00375016 

-0.00009449 

-0.00145509 

0.00028287 

G(10,  6) 

-0.00236521 

-0.00054373 

-0.00085616 

-0.00046020 

G(10,  7) 

-0.00159140 

0.00002960 

-0.00173427 

0.00052795 

G(10,  8) 

-0.00261024 

0.00041239 

0.00012128 

0.00057350 

G(10,  9) 

0.00000505 

0.00000053 

0.00000181 

-0.00000112 

G(10,10) 

-0.00000605 

0.00000012 

0.00000627 

-0.00000064 

Appendix  B 

Functional  Representations  G(n,m,a)  of  the  Expansion  Coefficients  of  the  Associated  Legendre 
Functions  Which  Represent  the  Magnetic  Potential  Due  to  Magnetopause  Surface  Currents 
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G(5,l,a)  =  -  (.5*(1.6*((.38/jc2)*Cos(a)  +  (.95/(97c2))*Cos(3a)  +  (.32/(257c2))*Cos(5a) 
(.4/(497t2))*Cos(7a))  -  .0985*Cos(a)  +  .004*Cos(4a)  +  .002*Cos(2a))  -  .005) 


G(5,2,a)  =  -  (.006*Cos(3a)  +  .008*Cos(2a)  -  .0034*Cos(a)  -  .002*Cos(4a)  -  .001*Cos(5a)  -  .0695) 
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G(7,2,a)  =  -  (1.03*(.0022*Cos(a)  -  .0075*Cos(2a)  +  .0034*Cos(3a)  -  .01*Cos(4a))  -  .036) 
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G(10,0,a)  =  -  (.7*(.00019»Cos(a)  -  .0004*Cos(2a)  +  ,0005*Cos(3a)  -  .00026*Cos(4a))  +  .000002) 


